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ABSTRACT 

Observations of the dust emission using millimetre/submillimetre bolometer arrays 
can be contaminated by molecular line flux, such as flux from 12 CO. As the brightest 
molecular line in the submillimetre, it is important to quantify the contribution of CO 
flux to the dust continuum bands. Conversion factors were used to convert molecular 
line integrated intensities to flux detected by bolometer arrays in mJy beam -1 . These 
factors were calculated for 12 CO line integrated intensities to the SCUBA-2 850 
and 450 |Ltm bands. The conversion factors were then applied to HARP 12 CO 3-2 
maps of NGC 1333 in the Perseus complex and NGC 2071 and NGC 2024 in the 
Orion B molecular cloud complex to quantify the respective 12 CO flux contribution to 
the 850 |Ltm dust continuum emission. Sources with high molecular line contamination 
were analysed in further detail for molecular outflows and heating by nearby stars to 
determine the cause of the 12 CO contribution. The majority of sources had a 12 CO 
3-2 flux contribution under 20 per cent. However, in regions of molecular outflows, 
the 12 CO can dominate the source dust continuum (up to 79 per cent contamination) 
with 12 CO fluxes reaching ~68 mJy beam -1 . 

Key words: submillimetre - instrumentation: detectors - dust, extinction - ISM: 
molecules - stars: formation - ISM: jets and outflows. 



1 INTRODUCTION 

Dust continuum is a useful tracer of star formation. The 
submillimetre (roughly 300 to 1000 \im) detection of the 
dust emission identifies Class and prestellar cores as well 
as disks, filamentary structure in molecular clouds, and the 
dust and gas masses of galaxies. In order to quantify the flux 
from dust in the submillimetre wavelengths, the heat gen- 
erated by the radiation is measured by bolometers with the 



detected wavelength range defined by wide-band filters ( Hoi 



land et al.||2002). These observations of broadband contin- 



uum emission from the dust can be contaminated by molec- 
ular line flux, particularly from 12 CO, which is the second 
most abundant molecule in the interstellar medium (after 
H2) with strong emission lines in the submillimetre ( |John 
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[ stone et a l. 2003; Seaq uist et al.|2004| [Hatchell Dunham| 



2009 



Zhu et al.||2003 ). Since the molecular line contamina- 



tion depends explicitly on the bandwidth and wavelength 
of the bolometer, it is important to quantify the potential 
contribution from molecular lines to make accurate flux mea- 
surements of the submillimetre dust emission used to calcu- 
late masses. 

The CO line contribution can be quantified by com- 
paring observations of the dust continuum emission and the 
CO line emis sion (|Gordon||1995|. Past research ( Johnstone 



&; Bally|1999l|Davis et al.|2000||Papadopoulos fe Allen| 2000 



Tothill et al. 2002) has studied the Submillimetre Com- 



mon User Bolometer Array (SCUBA) at the James Clerk 
Maxwell Telescope (JCMT), where line contribution from 
the 12 CO 3-2 line was found to range from little to tens 
of per cent in the 850 \im band. Other studies have exam- 
ined contamination in various bolometer instruments, in- 
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eluding MAMBO-II, Bolocam, and SHARC-II. The Submil- 
limetre High Angular Resolution Camera II (SHARC-II) op- 
erates at the same wavelength range as SCUBA (450 [im 
and 850 |^m), but also includes a 350 \im filter (780 to 
910 GHz). The 350 \un SHARC-II continuum could be po- 
tentially contaminated by the 12 CO 7-6 line (806 GHz) 
up to ^20 per cent, similar to the 12 CO 3-2 contamina- 



tion to the SCUBA 850 \im continuum ( |Hatchell Sz Purr 
ham|2009 ). While Bolocam (operated at Caltech at 1.1 mm 



with a 250 to 300 GHz filter) has been designed to ex- 
clude 12 CO line contamination, the Max Planck Millime- 
tre Bolometer II (MAMBO-II; operated by the Max Planck 
Institut fur Radioastronomie at 1.2 mm with a ^210 to 
290 GHz filter) includes 12 CO 2-1 (230 GHz) molecular line 
emission which could potentially increase flux at most a few 
percent ( Aguirre et al.|2011 ). Other possible contamination 
for Bolocam and MAMBO-II could result from other known 
molecular lines in clouds, including SiO 6-5 (^260 GHz) and 
HCN 3-2 (-258 GHz). 

The successor to SCUBA is SCUBA-2, a 10,000 pixel 
submillimetre camera on the JCMT with eight sub- arrays. 
SCUBA-2 operates at 450 [im and 850 (xm, like SCUBA, 
and can be susceptible to significant molecular line con- 
tamination. In regards to 12 CO, both SCUBA-2 bandpass 
filters have a central transmission peak near a 12 CO line: 
the 850 \xm bandpass filter centre is at 347 GHz near the 
12 CO 3-2 line at 345.796 GHz and the 450 \un bandpass fil- 
ter centre is 664 GHz near the 12 CO 6-5 line at 691.473 GHz. 
The proximity of the 12 CO line frequencies to the centres of 
the transmission peaks makes significant CO contamination 
in SCUBA-2 maps likely. 

In this paper, we have calculated conversion factors 
used to convert maps of molecular line integrated intensity 
(K km s _1 ) to maps of molecular line flux (mjy beam -1 ) 
that contaminates the dust continuum emission. These con- 
version factors were calculated for 12 CO 3-2 contributions 
to the 850 \xm SCUBA-2 dust continuum emission and for 
12 CO 6-5 contributions to the 450 fxm SCUBA-2 contin- 
uum. Conversion factors were applied to HARP 12 CO maps 
of NGC 1333, a region in the Perseus molecular cloud com- 
plex, and NGC 2071 and NGC 2024, regions in the Orion B 
molecular cloud complex, to calculate the contamination di- 
rectly by measuring fluxes and masses from a list of sources. 
Once the 12 CO contamination to the source fluxes was calcu- 
lated, the sources with the highest contamination were anal- 
ysed in more detail to determine the cause of the molecular 
flux contribution, e.g. molecular outflows or hot molecular 
gas from nearby stars. 

This paper has six sections. In Section [2] we give details 
of the calculation of the conversion factors considering dif- 
ferent weather grades (Section 2.1). Section [3] presents the 
resulting conversion factors calculated for 1 CO 3-2 (con- 
tributing to the 850 |^m continuum band) and 12 CO 6-5 
(contributing to the 450 \im continuum band). Section [4] 
explains how the conversion factors were applied to HARP 
12 CO 3-2 maps to quantify the amount of contamination to 
the 850 \xm SCUBA-2 dust continuum maps of three regions, 
NGC 1333, NGC 2071, and NGC 2024. Section [XT] and [O] 
introduce flux and mass calculations for a list of sources 
in the regions, and Section |4.3| presents further analysis of 
sources with high 12 CO contamination. Section [H] discusses 
the results, including the effects of environment and loca- 



tion (molecular outflows and hot ambient gas from nearby 
stars) on the 12 CO source contamination. This section also 
includes estimates of the 12 CO 6-5 line to the 450 \im 
SCUBA-2 dust emission. Lastly, Section [6] summarises the 
conclusions drawn from this work. 



2 METHOD 

Molecular line emission is typically measured as an intensity 
or surface brightness in terms of the Ray leigh- Jeans (R-J) 
brightness temperature (in Kelvin), while the dust contin- 
uum fluxes are given in Janskys measured over the telescope 
beam area (Jy beam -1 ). In order to convert 12 CO line inten- 
sities to pseudo-continuum fluxes, the intensity of a molecu- 
lar line must be converted into the flux of the line using the 
following relation 



-/ 



I dQ » in, 



(1) 



where / is the intensity and Q is the telescope beam area. 
The intensity is measured as a main-beam brightness tem- 
perature Tmb in Kelvin and converted to intensity using 



2v' 



kT M B 



2k 



(2) 



where v is the frequency, A is the wavelength, and k is the 
Boltzmann constant. 

A narrow molecular line within a filter contributes flux 
over a smaller frequency range (vn ne ) than continuum emis- 
sion across the filter. To obtain the flux from the molecular 
line, the average intensity (I) must be calculated over the 
full filter band, i.e. 



(i) 



J I u (line) g u (\me) dv 
J g v dv 



(3) 



where /^(line) is the intensity of the molecular line, g v {Ymo) 
is the filter passband (transmission) at the frequency of the 
molecular line, and J g v dv is the integrated filter passband 
(transmission) across the full range of filter frequencies. Us- 
ing Equation [5] and the Doppler shift, Av/v = Av/c, Equa- 
tion [3] can be converted to Tmb- 



(I) 



J I u (v) g u (v) dv _ 2kv 3 s^(line) 

C 3 J g u dv 



j g v dv 



dv, (4) 



where J Tmb is the velocity integrated main-beam bright- 
ness temperature, or integrated intensity. Using these calcu- 
lations of intensity, it follows from Equation [I] that 

F u _ 2kv 3 3^ (line 



mjy beam 



f g v dv 



Tmb dv 



(5) 



A similar calculation was used in Seaquist et al. (2004). 



The main-beam brightness temperature Tmb is used 
rather than the antenna temperature TX for analysing small- 
scale emission as long as the continuum calibration accounts 
for the same beam efficiencies as the molecular line emis- 
sion. This is true for the SCUBA-2 and HARP pairing at 
the JCMT. The beam efficiencies are discussed further in 
Section [4] TJ is related to Tmb by the following: 

T4 — ?7mbTmb, (6) 
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where ?7mb is the main-beam efficiency factor. The efficiency 
factor that takes into account emission at larger scales is dis- 
cussed in Sect ion [2^1 At ~345 G Hz with HARP on JCMT, 
?7mb is 0.61 (|Buckle et al.|[2009). The telescope beam area, 



also discussed further in Section [2^2] is measured in steradi- 
ans (sterad) and obtained from the full-width-half-maximum 
(FWHM) of a Gaussian beam using Qb = 2na 2 where 
the FWHM B 2\/2 In 2a: 



= ln2 V " J \] 



, i , (7) 

sterad 4 In 2 V " ) Vl80 x 36007 v ; 

Using Equation [5] a molecular line conversion factor, 
C, can be calculated to convert molecular line maps, mea- 
sured in the velocity integrated main-beam temperature 
J Tmb dv (K km s _1 ), into maps of line flux (mJy beam -1 ) 
that contributes to the observed continuum emission, 



C 



mJy beam 1 per K km s 



j Tmb dv 
2kv 3 ^(line) 
c 3 J g u dv 



(8) 



where frequencies are measured in GHz and 1 Jy = 
10 26 W m 2 Hz" 1 = 10 23 erg s" 1 cm" 2 Hz" 1 . Note that the 
beam size is wavelength dependent, where oc A. The dif- 
ference in beam size between the 12 CO line and the SCUBA- 
2 filter is not taken into account. 



2.1 Line Conversion Factors 

To calculate the conversion factor C from Equation [9] the 
SCUBA-2 filter profiles and the added atmospheric trans- 
mission were used to find j g v dv and ^(line). The SCUBA- 
2 850 |im and 450 |im filter profiles are shown in the bottom 
plot of Figures [l] and [2] The SCUBA-2 filter profiles are a re- 
sult of stacking all of the filters (thermal and bandpass filters 
as well as the cryostat window and dichroic) that form the 
continuum bandpasses when combined with the atmosphere. 
The bandpasses are the main filters defining the transmis- 
sion window, where passbands are the range of frequencies 
with a signal passing through the filter and stopbands define 
frequency ranges with a signal attenuated by the filter. The 
main infrared (IR; thermal) blocking filters are designed to 
block transmission at higher frequencies (IR and optical). 
For this study, a constant value for these filters has been as- 
sumed due to the high transmission in our frequency range. 
For further information, see the JCMT website regarding 
the cryostat window, filter and dichroic specification and 
measurement^ 

The JCMT has a system that describes the atmospheric 
conditions ranging from weather grades 1-5. The atmo- 
spheric conditions are based on precipitable water vapour 
(PWV) levels (in mm) that correspond to different sky opac- 
ities at 225 GHz, or T225- The relation between PWV and 
T225 is the following (JCMT Telescope Overview website^] 

T225 « 0.01 + (0.04 x PWV) (9) 

The JCMT weather grades are defined as: 



1 |http : //www . j ach . hawaii . e du/ JCMT/ continuum/ scub a2/ | 



• Grade 1: PWV: < 1 mm, T225: < 0.05 

• Grade 2: PWV: 1 to 1.75 mm, T225: 0.05 to 0.08 

• Grade 3: PWV: 1.75 to 2.75 mm, T225: 0.08 to 0.12 

• Grade 4: PWV: 2.75 to 4.75 mm, T225: 0.12 to 0.20 

• Grade 5: PWV: > 4.75 mm, T225: > 0.20 

The continuum bandpass transmission of the filter pro- 
file and atmosphere combined varies depending on atmo- 
spheric conditions. Therefore, the CO contamination was 
calculated based on these five weather grades. Plots of at- 
mospheric transmission corresponding to these conditions 
can be found in the top half of Figures [I] and [2] and are la- 
belled according to the respective water vapour levels (for 
more detail, see CSO Atmospheric Transmission Interactive 
Plotter website) The SCUBA-2 850 fxm and 450 fxm fil- 
ter profiles were multiplied by each individual atmospheric 
transmission profile to produce continuum bandpass profiles 
at each weather grade, shown in the bottom of Figures [I] and 

m 

In Equation 5 J g u dv is the integrated SCUBA-2 con- 
tinuum bandpass and is calculated as the the sum of g v x Sv 
(where Sv is 0.01 GHz) at each corresponding frequency 
with units in GHz. The transmission of 12 CO, ^(line), is 
the transmission of the SCUBA-2 850 \im continuum band- 
pass at 345.7960 GHz (the rest frequency of the 12 CO 3- 
2 line) and the transmission of the SCUBA-2 450 \im con- 
tinuum bandpass at 691.4731 GHz (the rest frequency of the 
12 CO 6-5 line). To calculate conversion factors for redshifted 
lines, the frequency v and transmission g u (\me) (from Fig- 
ures [l] & [2]) in Equation [9] must be changed appropriately. 
The SCUBA-2 beam size is calculated using Equation [7] as- 
suming the main-beam FWHM 0b is 13.8" at 850 |^m and 
8.3" at 450 \im. The possibility of an associated secondary 
beam is discussed in Section [Z2l 



2.2 Telescope Beam Area 

The beam profile of the original SCUBA instrument di- 
verged from a single Gaussian and displayed a beam profile 
of two combined Gaussians: a primary beam roughly cor- 
responding to t he assumed FWHM and a secondary beam 
of 40" FWHM ( |Di Francesco et al.p008) . For the 450 \un 
maps, the primary beam had a 8.5 FWHM with a 0.90 
relative amplitude and the secondary beam had a 0.10 rel- 
ative amplitude. For the 850 \im maps, the primary beam 
had a 13.5" FWHM with a 0.96 relative amplitude and the 
secondary beam had a 0.04 relative amplitude. 

It is possible that the total beam of SCUBA-2 also in- 
cludes a broader secondary component. As explained in Sec- 
tion |2T] the telescope beam areas for this study were cal- 
culated assuming FWHMs for the 450 fxm (8.3") & 850 \im 
(13.8") SCUBA-2 beams. The primary beam is appropri- 
ate for studying the CO contamination in compact sources 
(small-scale emission), e.g. protostars and small outflows. 
However, if the CO emission is both bright and extended, 
then it may be necessary to include the secondary beam 
in the calculation. By fitting a two component Gaussian to 
coadded SCUBA-2 maps of Uranus, the 450 ]im primary 
FWHM is 8.7" (relative amplitude 0.83) and secondary 



filter/ 

^~ www . j ach . hawaii . edu/ JCMT/overview/tel_overview 
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Figure 1. Upper: Plots of the atmospheric transmission at 300 to 400 GHz, given 0.5 to 4.75 mm of precipitable water vapour. Lower: 
The upper line is the profile of the SCUBA-2 850 urn filter, and the lines beneath represent the SCUBA-2 filter with the addition of 
the atmospheric transmission at varying water vapour levels. The 12 CO 3-2 line is plotted at 345.7960 GHz. As shown in Table [l] the 
atmospheric transmission corresponds to different bands of weather used for observations. 
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Figure 2. Upper: Plots of the atmospheric transmission at 550 to 750 GHz, given 0.5 to 4.75 mm of precipitable water vapour. Lower: 
The upper line is the profile of the SCUBA-2 450 \xm filter, and the lower lines represent the SCUBA-2 filter with the addition of the 
atmospheric transmission at varying water vapour levels. The 12 CO 6-5 line is plotted at 691.4731 GHz. 
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691.4731 
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Table 1. Line contribution factors for 12 CO lines in the SCUBA-2 850 |J.m Sz 450 u.m continuum bands. 



FWHM is 20.4" (relative amplitude 0.17) and the 850 [im 
primary FWHM is 13.9" (relative amplitude 0.97) and sec- 
ondary beam FWHM is 39.1" (relative amplitude 0.03) Q 
The effective FWHM becomes 11.6" and 15.3" for 450 urn 
and 850 urn beams respectively. This would cause the total 
beam area for 450 \im to be higher by a factor of 2.0 and 
the 850 |^m total beam area to be higher by a factor or 1.2. 

If it is necessary to incorporate the secondary beam into 
the beam area calculation, then a new conversion factor can 
be calculated using Equation [9] The conversion factors are 
directly proportional to the telescope beam area. Assuming 
C' is the the conversion factor with the inclusion of both 
a primary and secondary beam (Q' B ) and C is the relation 
shown in Equation [9] then it follows from Equation [9] 



C = C 



Mb _ r 



/ FWHM 7 
I FWHM 



(10) 



where FWHM 7 is the effective FWHM. For large-scale 
and extended emission on scales significantly larger than 
the size of the telescope beam (greater than 13.8", the 
850 urn FWHM), it is also more appropriate to use ?7f ss , 
the forward spillover and scattering efficiency, to calibrate 
the CO emission instead of the main-beam efficiency ?7mb. 
The forward spillover and scattering efficiency measures the 
amount of coupling to an extended source up to a 30' diame- 



ter (measured by observing the Moon). Equation 10 becomes 



C = C 



( FWHM 7 



V FWHM J rj {i 



tjmb 



(11) 



3 RESULTS 

The 12 CO conversion factors, C, for SCUBA-2 are 
listed in Table [l] The 12 CO 3-2 conversion factors (in 
mjy beam -1 per K km s -1 ) range from 0.63 (Grade 1) to 
0.77 (Grade 5) with a mid value of 0.70 (Grade 3). The con- 
version factors change depending on the atmospheric con- 
ditions that affect the continuum bandpass profile. Since 
each increase in T225 causes the 850 \im continuum bandpass 
profile to become narrower with less overall transmission, 
j g v dv in Equation [S] shrinks faster than the transmission of 
12 CO, ^(line). Therefore, the 12 CO 3-2 line contribution to 
the flux is lowest in Grade 1 weather and steadily increases 
with each step to Grade 5 weather. 

For the 450 urn continuum bandpass profile, the oppo- 
site trend is seen. The 12 CO line contribution to the 450 |^m 
flux is highest in Grade 1 but steadily decreases with each 
step to Grade 5 weather. In most cases, observations us- 
ing SCUBA-2 450 \xm would only be taken in Grade 1 
to 3 weather due to the decreased transmission in higher 
weather grades. The 12 CO 6-5 conversion factors range from 
0.64 (Grade 1) to 0.35 (Grade 5) with a mid value of 0.51 
(Grade 3). 

Note that the contamination is expected to have differ- 
ent behaviour between the 450 \im and 850 \im filters. The 
12 CO 3-2 line is in the centre of the 850 urn filter with the 
bulk of the transmission while the 12 CO 6-5 line is close to 
the edge of the 450 urn filter with lower transmission. The 
molecular line contribution to the 450 urn band decreases 
with weather grade because of the increasing attenuation of 
the 12 CO 6-5 line. 



The increase of telescope beam area caused from the inclu- 
sion of the secondary beam is somewhat counterbalanced 
by the use of = T^/r]{ ss rather than Tmb, accounting for 
the more efficient telescope coupling to large-scale emission 
(?7fss = 0.71 compared to ?7mb = 0.61). 

The secondary beam and potential changes in conver- 
sion factors for large-scale emission are further discussed in 
Section 



4 Measured primary beam sizes are slightly larger than the sizes 
quoted in Section |2.1| This is potentially due to small pointing 
shifts between coadded maps. 



4 APPLICATIONS TO OBSERVATIONS 

The conversion factors calculated in Section [3] were ap- 
plied to HARP 12 CO 3-2 maps and compared to SCUBA- 
2 850 um dust emission maps to measure the 12 CO contam- 
ination directly. Three different regions were used for this 
study: NGC 1333, NGC 2071, and NGC 2024. By quantify- 
ing the percentage of contamination to the dust continuum 
flux, we can determine regions more likely to be contami- 
nated by CO (i.e. regions with molecular outflows or nearby 
stars) . 
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4.1 Flux Calculations 

The continuum observations were taken with SCUBA-2 at 
450 pii and 850 \im in each region in 2010 during the 
SCUBA-2 Shared Risk Observing (S2SRO) campaign when 
SCUBA-2 had two science grade arrays (one at 450 \xm and 
one at 850 \im) installed. Observations were taken in Grade 
2 weather conditions. 

The Heterodyne Array Receiver Programme (HARP) is 
a 16 pixel array receiver (16 receptors separated by 30" and a 
footprint of 2'). The 12 CO 3-2 data cube for NGC 1333 was 
observed in J anuary 2007 u sing raster mapping techniques 
(for details see |Curtis et al.|2010| ). The 12 CO 3-2 data cubes 
for NGC 2071 and NGC 2024 were observed in November 



2007 using raster mapping techniques as well (see Buckle 



et 



"aLl |2010 > 



Both datacubes were rebinned to 0.42 km s~ 
velocity channels and converted to Tmb using a main-beam 
efficiency ?7mb of 0.61. 

Similar to other ground-based bolometer arrays, the 
limited, single-subarray version of SCUBA-2 available for 
S2SRO reproduced maps that are not sensitive to large-scale 
emission, in this case on scales larger than the single sub- 
array fleld-of-view (~4'). However, the HARP 12 CO maps 
still contain this large-scale flux. A simple application of 
the CO conversion factors from Table ^ to the HARP maps 
would retain the large-scale structure and overestimate the 
CO contamination. In order to account for the spatial fil- 
tering inherent in bolometer array reconstruction in a sim- 
ple way and subtract the large-scale flux from the maps, a 
Gaussian smoothing mask (GSM) filter was applied to both 
the 12 CO HARP integrated intensity maps (in K km s _1 ) 
and the 850 \im maps (mJy beam -1 ) of the regions. The 
GSM filter was designed to minimise emission from struc- 
ture on scales inaccessible to SCUBA-2 at this time. For 
the HARP maps, GSM filters were created by convolving a 
HARP 12 CO contamination map directly with a Gaussian a 
few arc minutes FWHM in size and subtracting the result- 
ing smoothed map from the original map. For the SCUBA- 
2 maps, it was necessary to first create a thresholded map 
for masking bright protostars and convolve the thresholded 
map with a Gaussian the same FWHM in size. The resulting 
smoothed map was then subtracted from the original map. 

A l' FWHM Gaussian was chosen for generating GSM 
maps after analysing l' to 3 ; GSM filter sizes, further dis- 
cussed in Section |4.1.4| Figure [3] shows examples of the 
SCUBA-2 GSM processed maps for protostellar cores LBS- 
MM18 (NGC 2071-IRS) in NGC 2071 (see 



2001), FIR 1-7 in NGC 2024(see Richer et al 



et al ||2010| , and SVS13 in NGC 1333 (see |Hatchell et al. 
2007). The brightest CO features in the maps are the re- 



Motte et al. 



1989 Buckle 



suit of molecular outflows driven by the dense cores in the 
regions. Dust emission contours in the outflows of these re- 
gions clearly follow the 12 CO 3-2 emission, indicating the 
CO contamination is strong enough to be directly detected 
in the dust continuum. 

To study the CO contamination quantitatively, aper- 
ture photometry with a 15" radius was applied to lists of 
known submillimetre sources, listed in Tableland further 



a sky background of zero with flux uncertainties based on 
the sky RMS and include a correction for the Gaussian beam 
(Eno ch et al.|2006 ). Therefore, a point-source has the same 
integrated flux density in any size aperture. 



4.1.1 Application to NGC 1333 

NGC 1333 is a reflection nebula in the Perseus molecular 
cloud and is characterised by early stage star formation 
of age less than 1 Myr ( |Lada et al.||1996| |Wilking et al. 
|2004| ). The flux calibration for the S2SRO maps of this re- 
gion was the CRL618 nebula and pointing checks were from 
the active galactic nucleus 3C84. A flux conversion factor 
(FCF) of 500 Jy beam -1 pW" 1 was used for NGC 1333 to 
convert the maps into mJy beam -1 . Sources were chosen 
from a list of cores in NGC 1333 ( |Hatchell et al.|2007| that 
had been previously identified in the submillimetre using 
SCUBA (iHatchell et al.||2005| and Bolocam ([Enoch et al 



2006 ) with a total of 35 sources in the area covered by the 



SCUBA-2 map. These sources include a mixture of proto- 
stellar and starless cores. For further information regarding 



HARP observations, see Curtis et al. (2010). 



Figure [4] shows the source fluxes from the SCUBA-2 
850 \xm and 12 CO 3-2 Grade 2 contamination maps and 
the percentage contribution of 12 CO 3-2 flux to 850 \im 
SCUBA-2 flux. It should be noted that all of the sources 
have 12 CO contributions less than 20 per cent and every 
source except one (source 21) has a contribution less than 
10 per cent. Source 21 is further discussed in Section |4~3[ 



4.1.2 Application to NGC 2071 

NGC 2071 is a region in the Orion B molecular cloud. 
Sources were chosen from a list of young stellar objects 



that had been previously identified using SCUBA (Nutter 



&; Ward-Thompson||2007| with a total of 50 sources in the 



area covered by the SCUBA-2 map. A flux conversion fac- 
tor of 685 Jy beam -1 pW" 1 was used for NGC 2071 as well 
as NGC 2024 (FCF value valid for October 2010 reduction, 
equivalent to Nutter et al, in prep.). For further informa- 
tionon HARP observations of NGC 2071, see |Buckle et al.| 
( [2010 ). 

Figure [5] shows the source fluxes from the SCUBA-2 
850 \im and 12 CO 3-2 Grade 2 contamination maps and 
the percentage contribution of 12 CO 3-2 flux to SCUBA-2 
850 |^m flux. Note that the majority of sources have 12 CO 
contributions of less than 20 per cent in Grade 2 weather. 
Four sources (sources 70, 74, 88, and 91) have 12 CO con- 
tributions greater than 20 per cent (ranging from 34 to 
79 per cent). Sources with a higher 12 CO contamination are 
further discussed in Section f4. 31 



4.1.3 Application to NGC 2024 



4.1.1 



4.1.2 



discussed in Sections 
radius was chosen based on source proximity and the possi- 
bility of the aperture diameter extending to a neighbouring 
source. Integrated flux densities are calculated by assuming 



NGC 2024 is another emission nebula in the Orion B molec- 
ular cloud. Sources were chosen from a list of young stellar 
objects in NGC 2024 that had been previously observed us- 
and |4.1.3| A15" aperture ing SCUBA ( [Nutter fe Ward-Thompson||2007| ) with a total 



of 24 sources in t he are a covered by the SCUBA-2 map. As 
stated in Section 4.1.2 a FCF of 685 Jy beam -1 pW -1 was 



used to correspond with current studies of Orion B (FCF 
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RA (J2000) 

Figure 3. Regions in the SCUBA-2 850 urn GSM-filtered maps of NGC 2071, NGC 2024, and NGC 1333 where 12 CO emission contributes 
strongly to the 850 urn flux. Black contours correspond to the SCUBA-2 850 u.m dust continuum maps and white contours correspond 
to HARP 12 CO 3-2 contamination maps. Top Left: Close-up of LBS-MM18 (NGC2071-IRS) and corresponding outflow. Both sets of 
contours correspond to flux at 20 and 45 mJy beam -1 . Noticeable 12 CO flux contribution in the lower right lobe of the molecular outflow. 
Top Right: Close-up of FIR 1-7 and corresponding outflow. Both contours correspond to flux at 10, 20, and 45 mJy beam -1 . Noticeable 
12 CO flux contribution in the lower lobe of the molecular outflow. Bottom: Close-up of SVS13 and corresponding outflow. Both sets of 
contours correspond to flux at 10, 20, and 35 mJy beam -1 . Noticeable 12 CO flux contribution in the right lobe of the molecular outflow. 



value valid for October 2010 reduction, equivalent to Nutter 
et al, in prep.). For further information on HARP observa- 



tions of NGC 2024, see Buckle et al. 



(2010). 



Figure [6] shows the source fluxes from the SCUBA-2 
850 \im and 12 CO 3-2 Grade 2 contamination maps and 
the percentage contribution of 12 CO 3-2 flux to SCUBA-2 
850 |^m flux. Note that the majority of sources have 12 CO 
contributions of less than 20 per cent. One source (source 48) 
has a 

contamination 



12 CO contribution more than 20 per cent (34 per cent 
further discussed in Section 



4.3 



4.1.4 GSM Analysis 

For SCUBA-2, GSM filters were created by applying upper 
thresholds to the original SCUBA-2 maps which acted as 
a mask for source emission and convolving the thresholded 
maps with a Gaussian a few arcminutes in FWHM size (sim- 
ilar to |Reid fe Wilson|2005| |Kirk et al.|2006| ); this is a stan- 
dard technique for the SCUBA-2 data. Negative regions of 
flux, known as negative 'bowls,' surround very strong sources 
in the SCUBA-2 maps and are produced in the map re- 



construction process (see Johnstone et al.||2000| . SCUBA-2 
thresholding was necessary to prevent introducing new nega- 
tive bowls in the image caused by smoothing and subtracting 
bright continuum sources. Without thresholding, artificial 
negative bowling would have been further added to the map, 
causing negative flux to lower source fluxes and increase the 
calculated 12 CO 3-2 contamination. To test the effects of 
changing the filter size, S2SRO maps were analysed with 
an upper threshold of 15 mJy beam -1 for NGC 1333 and 
30 mJy beam" 1 for NGC 2071 and NGC 2024 (used to mask 
out bright sources) and l' to 3' FWHM Gaussian smooth- 
ing. Aperture fluxes from each of the GSM maps (l x , 2', and 
3') and the original, unsmoothed S2SRO maps were found 
to agree within uncertainties, indicating that the emission 
on scales which would have been affected by the GSM filter 
had been filtered out by the SCUBA-2 map reconstruction. 
The l' and 2' GSM filters were further analysed with ap- 
plication to the HARP 12 CO maps based on the similarity 
between S2SRO and SCUBA maps, for which scales greater 
than 2' are known to be poorly reproduced ( Hatchell et al.| 



2007). 
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Figure 4. Top: The SCUBA-2 source fluxes calculated from the 850 u.m continuum and 12 CO 3-2 contamination maps (Grade 2 weather) 
of NGC 1333. Note several sources (26, 42, and 66) have 12 CO flux contributions of mjy beam -1 . Bottom: The percentage contribution 
to the SCUBA-2 fluxes from the 12 CO contamination maps. Numbers are given arbitrarily to the sources and were based on the original 
list of SCUBA and Bolocam cores ( |Hatchell et al.|2007| >. 



For the HARP 12 CO maps, thresholding was not re- 
quired because 12 CO mainly traces molecular outflows which 
have bright, extended structures on scales not fully repro- 
duced by SCUBA-2. To generate the 12 CO GSM maps, GSM 
filters with 1' and 2' FWHM Gaussians were directly applied 
and subtracted from the original maps. Negative flux regions 
in the final HARP maps resulting from the oversubtraction 
of background flux estimated by the GSM filter was set to 
mjy beam -1 to prevent biasing the source fluxes in the 
aperture photometry process. The S2SRO V and 2' GSM 
maps were subtracted from the corresponding 12 CO GSM 
maps (Grade 2 contamination) to analyse the GSM filter ef- 
fectiveness in matching the spatial filtering of the SCUBA-2 
maps. Positive flux in the residuals indicates higher 12 CO 
flux than 850 \im dust continuum flux, suggesting that the 
GSM filter size should be scaled down to subtract smaller 
scale emission. In each of the regions, the 2! GSM residuals 
were found to overestimate the 12 CO flux contribution to 
the dust continuum. On average, the 2' GSM map residuals 
were 1.4 to 1.5 times greater than the 1' GSM. The 1' GSM 
filters were applied to both the S2SRO and HARP maps for 
consistency in eliminating flux on scales of l' and above. 



With the full complement of subarrays, SCUBA-2 is 
likely to recover more large-scale structure and continuum 
fluxes may increase further. A comparison of 12 CO contam- 
ination on scales of l' or greater will have to wait for full 
SCUBA-2 operations. 



4.2 Mass Calculations 

The calculation of the dust continuum flux from pre- and 
protostellar sources in a molecular cloud can be used to ob- 
tain source masses (Hildebrand 1983| ). Depending on the 
molecular cloud environment surrounding the sources, con- 
tamination from the 12 CO line emission may affect low- and 
high-mass sources, leading to a varying level of source con- 
tamination. Therefore, the masses of sources were calculated 
using the relation between the dust and gas mass and the 
total source dust continuum flux (e.g. |Seaquist et al.| |2004 
Enoch et al.|2006| ), 



M 



S850D 2 



^850^850 (Td) ' 



(12) 



where Ssso is the flux from 15" radius aperture photometry 
at 850 |^m, D is the distance to the source, ftsso is the dust 
opacity at 850 |^m and Bs5o(Td) is the Planck function at 
850 \im for the dust temperature Td. 



4.2.1 Mass Calculations for NGC 1333 

For NGC 1333, a distance of 250 pc was assumed for the 
mass calculations. Hatchell et al. (2007) used a distance of 

factor of 1.6. A 



320 pc, which would increase masses by 
temperature of 10 K was used as an estimate of the dust 
temperature, where dense regions that do not have internal 
heating are colder on the inside and warmer on the outside 
(Eva ns et al.|2 001). Cores with internal heating are warmer 
in the inner regions. For example, Class and Class I pro- 
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Figure 5. Top: The SCUBA-2 source fluxes calculated from the 850 u.m continuum and 12 CO 3-2 contamination maps (Grade 2 weather) 
of NGC 2071. Note source 67 has a 12 CO flux contribution of mjy beam -1 . Bottom: The percentage contribution to the SCUBA-2 
fluxes from the 12 CO contamination maps. Numbers are given arbitrarily to the sources and were based on the original list of cores from 
|Nutter fe Ward-Thompson| ( |2007| >. 
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Figure 6. Top: The SCUBA-2 source fluxes calculated from the 850 p.m continuum and 12 CO 3-2 contamination maps (Grade 2 weather) 
of NGC 2024. Bottom: The percentage contribution to the SCUBA-2 fluxes from the 12 CO contamination maps. Numbers are given 
arbitrarily to the sources and were based on the original list of cores from |Nutter fe Ward-Thompson| ( |2007| >. 
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tostars are found from models to have ~ Td = 15 K ( Shirley 
et al.|[2002 Young et al.||2003 ), but most of the dust mass 
is found in areas of lower temperatures. To cover pre- and 
protostellar sources, Td = 10 K is a commonly used aver- 
age. It should be noted that this value can overestimate the 
masses of protostellar sources by a factor of 2 to 3 if the 



temperature is warmer ( Enoch et al.|2006 ). 

The dust opacity ftsso is also uncertain for individual 
regions. Hatchell et al. ( 2007| assumed a dust opacity of 
0.012 cm g~ for the 850 SCUBA dust emission maps 



of NGC 1333 based on a gas/dust ratio of 161 (see Ossenkopf 
&; Henning|1994| . This dust opacity is at the low end of the 



assumed values and a dust opacity of 0.02 cm 2 g -1 at 850 \im 
could have been used (Kirk et al. 2006). Here, we choose a 
dust opacity of 0.012 cm g~ . If 0.02 cm 2 g -1 were used, 
then our masses would decrease by a factor of 1.7. 

Masses were calculated from source fluxes obtained 
from continuum emission with and without 12 CO contami- 
nation taken into account. The bottom portion of Figure [7] 
shows the ratio between these masses. Due to the CO flux 
contamination, the calculated source masses are being over- 
estimated by up to a factor of 1.2. 

4.2.2 Mass Calculations for NGC 2071 and NGC 2024 

For NGC 2071 and NGC 2024, parameters from past mass 
estimates ( |Nutter &; Ward- Thompson] |2007| were used to 
calculate source masses. A distance of 400 pc was assumed 
for both regions ( Brown et al.| [l994) and a temperature of 
20 K was assumed as an estimate of the dust temperature 
dLaunhardt et al.||1996j |Mitchell et~aLl|2001| |Johnstone k\ 



Bally 2006). An 850 \xm dust opacity of 0.01 cm g was 



used ((Andre et al.|2003||Andre et al.|1996|[Ward-Thompson 



et al. 11999 ), similar to the dust opacity used for NGC 1333. 



As in NGC 1333, the assumption of a single temperature for 
each source does introduce a potential bias in the masses. 
If 10 K was assumed, as for NGC 1333, then masses for 
NGC 2071 and NGC 2024 sources would be larger by a factor 
of 2. 

Masses were calculated from the source fluxes of the 
850 |^m continuum emission (see Section [4| with and with- 
out 12 CO 3-2 contribution taken into account. The bottom 
portion of Figures [8] and [5] show the ratio between these 
masses for NGC 2071 and NGC 2024 respectively. Due to 
the CO flux contamination, the calculated source masses are 
being overestimated by a factor up to 4.8 for NGC 2071 and 
1.5 for NGC 2024. 

4.3 Molecular Outflow Analysis 

The location of protostellar sources can help identify the 
potential causes of 12 CO contamination. The presence of 
protostellar molecular outflows and hot ambient gas from 
nearby stars results in bright 12 CO emission, making regions 
with these characteristics rife with CO contamination. 

Sources with high 12 CO contamination were examined 
in further detail using the HARP data cubes. In each of the 
three regions, sources with greater than 20 per cent con- 
tamination were defined as 'sources with high contamina- 
tion.' For NGC 1333, there were no sources with greater 
than 20 per cent contamination, excluding it from this por- 
tion of the high contamination study. In NGC 2071, four 



sources fulfilled the high contamination criterion and one 
source in NGC 2024 fulfilled the criterion. In order to iden- 
tify the cause of high contamination, the 12 CO spectra were 
extracted and analysed for molecular outflows. 

For NGC 2071 and NGC 2024, the linewing criterion 
used to identify a molecular outflow candidate was a linew- 
ing above 1.5 K (TJ) at ± 4 km s _1 from the core velocity, 
^lsr- This linewing c riterion method follows the method in 



Hatchell et al. 



(2007|. A core velocity of 10 km s 1 was used 
for all of the so urces in NGC 2024 a nd NGC 2071 based on 
C 18 3-2 data ( [Buckle et aT|2010 ). Linewing criteria were 
based on T% RMS values for the regions (5a). 

The linewing criterion identifies not only protostars 
driving molecular outflows, but also sources which are con- 
taminated by outflows along the line of sight. Outflow candi- 
dates were identified using the above criteria and examined 
further to determine if the source or another protostar was 
the outflow driving source. Sources with high contamination 
that were not outflow candidates were further analysed to 
determine if there were other causes behind the 12 CO con- 
tamination, such as a nearby star heating the gas. 

Highly contaminated sources are listed in Table [2] Ta- 
ble [2] includes the region, source number, RA and Dec, flux 
calculated from aperture photometry in mjy, 12 CO flux con- 
tamination in Grade 2 weather in mjy, percentage contam- 
ination from 12 CO, core velocity vlsr in km s _1 , and the 
final molecular outflow candidate result. The outflow nam- 



ing convention follows Hatchell Sz Dunham (2009), where a 
'y' is given when an outflow is present and 'n' is given when 
an outflow is not present. Sources are marked '?' when there 
is confusion as to the source of the outflow. In this case, the 
potential source causing the outflow detection is listed in a 
footnote. 



4.3.1 Sources in NGC 1 333 

Figure [To] shows the 850 \im SCUBA-2 map with blue con- 
tours tracing the blueshifted 12 CO HARP intensity J T% dv 
(integrated from -2.5 to 4.5 km s _1 ) and red contours trac- 
ing the redshifted 12 CO intensity (integrated from 10.5 to 
16.5 km s -1 ). Sources in NGC 1333 are denoted by the per- 
centage contamination, where 'x' denotes sources with 
to 10 per cent contamination and '+' denotes sources with 
10 to 20 per cent contamination. Source 21 had the high- 
est percentage contamination at 12 per cent. According to 
the linewing criteria used to identify a molecular outflow 
candidate (linewing above 1.5 K for TJ at ± 3 km s _1 
from the core velocity 7.9 km s _1 , following the criterion 
for NGC 1333 used in [Hatchell et al.||2007) , source 21 is a 
molecular outflow candidate that could potentially be the re- 
sult of a source at (J2000) 03:29:03.2, 31:15:59.0 (SVS13) or 
source at ( J200 0) 03:29:08.8, 31:15:18.1 (SK-16) ( |Hatchell & 
bunham|2009 ). A 12 CO spectrum for this source is included 
in Figure [Al] 



4.3.2 Sources in NGC 2071 

Figure [TT| shows 850 |^m SCUBA-2 map with blue contours 
tracing the blueshifted 12 CO HARP intensity j T\ dv (in- 
tegrated from -2.0 to 6.0 km s _1 ) and red contours trac- 
ing the redshifted 12 CO intensity (integrated from 14.0 to 
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Figure 7. Top: Mass calculations (in solar masses) of the different sources in NGC 1333. The masses were calculated using the SCUBA-2 
850 \im map and then recalculated excluding the flux contribution from 12 CO in different atmospheric conditions. Bottom: Ratio of the 
masses calculated from the flux without to with the 12 CO contribution. In both plots, uncertainties are calculated only from the source 
fluxes and do not include absolute calibration uncertainties. 



Table 2. List of sources categorised with high 12 CO contamination to the 850 (xm dust continuum. 



Region 


Source 


RA 
(J2000) 


Dec 
(J2000) 


Flmja] 

mJy 


12 CO (Grade 2) 
mJy 


Percentage Cont. 
(Grade 2) 


^LSR 
km s _1 


Outflow? 


NGC2071 


70 


05:46:57.6 


00:20:09 


495 ± 40 


224 ± 6 


46 ± 5 


10.0 


y 




74 


05:47:01.0 


00:20:42 


268 ± 40 


196 ± 6 


73 ± 13 


10.0 






88 


05:47:06.7 


00:23:14 


112 ± 39 


38 ± 3 


34 ± 14 


10.0 






91 


05:47:08.9 


00:23:56 


166 ± 40 


132 ± 5 


79 ± 22 


10.0 




NGC2024 


48 


05:41:19.9 


-01:54:16 


168 ± 33 


56 ± 3 


34 ± 9 


10.0 


y 



a SCUBA-2 fluxes appear to be lower than seen by SCUBA due to the subtraction of large-scale flux by SCUBA-2. 

b Could be due to a large, central blue outflow from source at (J2000) 05:47:06.9, 00:22:39 (source 84; LBS-MM19) or a sour ce at (J2000) 
05:47: 04.1, 00:21:58 (LBS-MM18; NGC2071-IRS), where LBS-MM18 was found to be responsible for driving the outflow in |Motte etaT] 
(2001) and both sources are confirmed Class 1 protostars detected using IRAC ( Nutter fc Ward-Thompson|2007| > . 



22.0 km s" 1 ). Sources in NGC 2071 are denoted by the 
percentage contamination, where 'x' denotes sources with 
to 10 per cent contamination, '+' denotes sources with 
10 to 20 per cent contamination, and 'O' denotes sources 
with greater than 20 per cent contamination. According to 
the linewing criteria used, all four sources with high con- 
tamination are molecular outflow candidates. Even though 
sources 70 and 74 have clear blue- and redshifted spectral 
linewings, all four sources appear to trace a large central out- 
flow that could be driven by a source at (J2000) 05:47:06.9, 
00:22:39 (source 84; LBS-MM19) or a source at (J2000) 
05:47:04.1, 00:21:58 (source 79; LBS-MM18; NGC2071-IRS), 
where LBS-MM18 was found to be responsible for driving 
the outflow in Mot te et al.| ( |2001| and both sources are con- 



firmed Class 1 protostars detected using IRAC (Nutter & 
Ward-Thompson||2007t . Note that sources 84 (LBS-MM19) 
and 90, both with 18 per cent contamination, also correlate 
with the central outflow. The 12 CO spectra for these sources 
are displayed in Figure [Al] 



4.3.3 Sources in NGC 2024 

Figure [l2| shows 850 \im SCUBA-2 map with blue contours 
tracing the blueshifted 12 CO HARP intensity J T\ dv (in- 
tegrated from -2.0 to 6.0 km s _1 ) and red contours trac- 
ing the redshifted 12 CO intensity (integrated from 14.0 to 
22.0 km s -1 ). Sources in NGC 2024 are denoted by the per- 
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Figure 8. Top: Mass calculations (in solar masses) of the different sources in NGC 2071. The masses were calculated using the SCUBA-2 
850 u.m map and then recalculated excluding the flux contribution from 12 CO in different atmospheric conditions. Bottom: Ratio of the 
masses calculated from the flux without to with the 12 CO contribution. In both plots, uncertainties are calculated only from the source 
fluxes and do not include absolute calibration uncertainties. 



centage contamination, where ' x ' denotes sources with to 
10 per cent contamination, '+' denotes sources with 10 to 
20 per cent contamination, and 'O' denotes sources with 
greater than 20 per cent contamination. According to the 
linewing criteria used, the single source with a high contami- 
nation is a molecular outflow candidate. The 12 CO spectrum 
for this source is listed in Figure [Al] 



5 DISCUSSION 

Typical 12 CO contamination levels in the observed SCUBA- 
2 850 \im emission from NGC 1333, NGC 2071, and 
NGC 2024 are under 20 per cent (this includes 95 per cent of 
sources, and 88 per cent of all sources have under 10 per cent 
contamination). Simila r results were found for SCUBA, 
(|2003|) suggested that 12 CO line con- 



where Johnstone et al 



tamination is typically under 10 per cent for submillimetre 
sources in Orion and |Davis et al. ( 2000| suggested contam- 
ination was ~10 per cent near the source V380 Orion NE. 

In locations where molecular outflows are present, 12 CO 
contamination can rise above 20 per cent and dominate 
the dust continuum (up to 79 per cent contamination), 
corresponding to a CO contribution ranging from 16 to 
68 mJy beam -1 for the sources analysed in this study. Peak 
12 CO fluxes found in the molecular outflows of NGC 1333, 
NGC 2071, and NGC 2024 maps reach even higher fluxes of 
84 mJy beam -1 , 154 mJy beam -1 , and 94 mJy beam -1 re- 
spectively. Our study suggests that molecular outflows can 



influence line contamination in sources with both high and 
low continuum fluxes (~ 100 to 500 mJy). This result agrees 
with that of Johnstone et al. (2003), where they concluded 



that the areas with warmer molecular gas temperatures and 
higher velocities (i.e. shocks and molecular outflows) were 
the only locations where 12 CO emission dominated the dust 
continuum flux due to the higher molecular line integrated 
intensities associated with such regions. In addition, Gueth 
|et al.| ( |2003| found a 20 per cent contamination for the well- 



known outflow source L1157. Hatchell &; Dunham (2009) 



similarly found a 20 to 30 per cent CO contamination level 
in the IRAS 03282+3035 outflow in Perseus. 

The 12 CO contamination combined with contamination 
from other molecular lines allows outflows to potentially be 
seen in continuum maps with a similar appearance as pro- 
tostellar cores or filamentary structure, which may be the 
case for SVS13 in NGC 1333 and the large, central out- 
flows in NGC 2071 and NGC 2024 that have strong evi- 
dence of molecular outflow lobes detected in the dust contin- 
uum emission maps (see Figure [3|. Other studies have sug- 
gested that regions involving molecular outflows can reach 
50 per cent 12 CO contamination, e.g. in the extended out- 
flow lobes of V380 Orion NE ( |Davis et aL|2000) , and even 
up to 100 per cent contamination, e.g. the central blue 
outflow region i n NGC 2071 drive n by source LBS-MM18 
(NGC2071-IRS; |Motte etHpSoT ). 



The FCF uncertainty from calibrator observations for 
the 850 pn S2SRO maps is 18 per cent (SMURF SCUBA-2 
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Figure 9. Top: Mass calculations (in solar masses) of the different sources in NGC 2024. The masses were calculated using the SCUBA-2 
850 fim map and then recalculated excluding the flux contribution from 12 CO in different atmospheric conditions. Bottom: Ratio of the 
masses calculated from the flux without to with the 12 CO contribution. In both plots, uncertainties are calculated only from the source 
fluxes and do not include absolute calibration uncertainties. 



Data Reduction Cookbook J] The calibration uncertainty of 
HARP observations at JCMT is estimated to be 20 per cent 
by |Buckle et al. ( 2009| . With contamination levels to 
SCUBA-2 850 ujxl less than or equal to 20 per cent for the 
majority of the sources, the typical contamination is less 
than or equivalent to the calibration uncertainty. Problems 
arise when the contamination is greater than than calibra- 
tion uncertainties, contributing a significant portion of flux 
and potentially dominating the dust continuum. For bright 
sources, it should be possible to use the subtracted 12 CO 
background to estimate the column density and hence the 
potential CO contamination to the 850 ujn SCUBA-2 dust 



emission, as suggested by Tothill et al. (2002). 



For faint sources, the insensitivity of SCUBA-2 to the 
large-scale dust emission introduces additional uncertain- 
ties. The spatially filtered maps created here do not appear 
to entirely subtract the total large-scale flux detected by 
HARP. Excess 12 CO flux seen as positive flux in residuals 



(Section 4.1.4) point to the need for a more detailed model 
of SCUBA-2 structure response. The large-scale reconstruc- 
tion issues limit analysing contamination where the SCUBA- 
2 flux is faint. This analysis is something we will consider 
for the full SCUBA-2 array, but not for the limited S2SRO 
data. If the SCUBA-2 and HARP maps were subtracted in 
order to account for the molecular line contamination, there 



5 |http : //star-www . rl . ac . uk/star/docs/scl9 . htx/no de40 . 
[html | 



is a possibility that an overcompensation for the molecular 
line flux would occur, creating regions of negative flux in the 
dust continuum map. 

An additional uncertainty in the contamination due to 
large-scale CO emission is the inclusion of the secondary 
beam in the calculation of the telescope beam area in Equa- 
tion [9] which increases the beam area by a factor of 1.2 
at 850 \im (determined from a new FWHM calculated in 
Section |2.2| ). Large-scale emission couples to the telescope 
beam as the efficiency factor ?7f ss instead of ?7mb (as in 
Equation 11). Therefore, the conversion factor for 12 CO 3- 



2 would increase by a factor of 1.1, which is insignificant 
compared to other uncertainties. The scales of large-scale 
emission that could cause significant signal without being 
taken out by the S2SRO common-mode subtraction range 
from 13.8" (the FWHM of the 850 [im beam) to l' (the 
Gaussian FWHM used in the GSM masking process). For 
the full SCUBA-2 array, CO emission on scales up to 8' in 
size could contribute. 

CO is not the only possible contributor in the 850 
and 450 \im bands. Studies of other molecular lines found 
SCUBA 850 nm line contamination from HCN, HNC, CN 
and methanol add together to form roughly 40 per cent of 
the total line contamination when observing other more en- 
ergetic sources, like the shocked region SK1-OMC3 ( |John-| 
stone et al.||2003 ) . Similar contamination was found in the 



Kleinmann-Low nebula, from SO and SO2 emission that was 
28 to 50 per cent of the total line contamination at 850 \im 
( |Serabyn fe Weisstein|[l995| |Groesbeck et al.||1994| ). Other 
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Figure 10. SCUBA-2 850 fim map of NGC 1333. The colour 
bar represents flux in mJy beam -1 . Blue contours correspond 
to blueshifted 12 CO 3-2 HARP intensity J dv (integrated 
from -2.5 to 4.5 km s -1 ). Red contours correspond to redshifted 
intensity (integrated from 10.5 to 16.5 km s -1 ). Contour levels are 
5, 10, 15, 25, 45, 65, and 85 K km s -1 . Sources in NGC 1333 are 
denoted by percentage contamination, where ' x ' denotes sources 
with to 10 per cent contamination and '+' denotes sources with 
10 to 20 per cent contamination. 
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Figure 12. SCUBA-2 850 \im map of NGC 2024. The colour 
bar represents flux in mJy beam -1 . Blue contours correspond 
to blueshifted 12 CO 3-2 HARP intensity J T\ dv (integrated 
from -2.0 to 6.0 km s -1 ). Red contours correspond to redshifted 
intensity (integrated from 14.0 to 22.0 km s — 1 ). Contour levels 
are 5, 10, 15, 25, 45, 65, 85, 105, 125, and 145 K km s -1 . Sources 
in NGC 2024 are denoted by percentage contamination, where ' x ' 
denotes sources with to 10 per cent contamination, denotes 
sources with 10 to 20 per cent contamination, and 'O' denotes 
sources with greater than 20 per cent contamination. 
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Figure 11. SCUBA-2 850 M-m map of NGC 2071. The colour 
bar represents flux in mJy beam -1 . Blue contours correspond 
to blueshifted 12 CO 3-2 HARP intensity J dv (integrated 
from -2.0 to 6.0 km s -1 ). Red contours correspond to redshifted 
intensity (integrated from 14.0 to 22.0 km s -1 ). Contour levels 
are 5, 10, 15, 25, 45, 65, 85, 105, 125, and 145 K km s -1 . Sources 
in NGC 2071 are denoted by percentage contamination, where ' x ' 
denotes sources with to 10 per cent contamination, denotes 
sources with 10 to 20 per cent contamination, and 'O' denotes 
sources with greater than 20 per cent contamination. 



studies have found the total line contamination by other 
molecular lines to be a factor of 2 to 3 times that from CO 
in outflows flTothill et aI1|2002| |Gueth et al.||2QQ3| ). Since 
molecular line contamination from other molecules is also 
likely, some features with low flux in the dust continuum 
may entirely be the result of line emission. 



5.1 



CO 6—5 Contamination 



For the 450 \im band, we have no 12 CO 6-5 maps with which 
to estimate the CO contamination directly. Using the line 
intensities from the 12 CO 3-2 HARP maps, we can instead 
predict the potential line contamination from 12 CO 6-5 to 
the SCUBA-2 450 |^m dust continuum signal. 

Assuming local thermodynamic equilibrium, we can 
estimate the ratio of the main-beam brightness temper- 
atures T M b for 12 CO 6-5 and 12 CO 3-2. We assume 
the excitation temperature, T ex , is equal to the kinetic 
temperature of the region, and is therefore the same for 
both 12 CO 3-2 and 12 CO 6-5. We also assume the par- 
tition function Z « 2T ex /To and the Gaussian line shape 
6>(peak) = 2c V2 In 2 / v Avv / 27r, yielding the relation (in 
CGS units): 



M 2 (J+1) 2 



~3h 

2V2 In 2 
Ava/2^ 



2T e , 



exp 



-(J+1)(J + 2)T 

2 Tex 



(13) 



where h is Planck's constant, \i is the permanent electric 
dipole moment of the molecule, J is the lower rotational level 
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of a linear molecule, and To is the gr oun d-state temperature 
(hvo/k) at 5.5 K. Using Equation 13 the ratio 12 CO 6- 
5/ 12 CO 3-2 is: 



T MB (6^5) 
T M b(3^2) 



(6) 2 exp ( 



-21 T > 

^exp(^) 
= 4exp(-15g-) 



(14) 



where J(6 5) = 5 for T MB (6 5) and J(3 2) = 2 for 
T MB (3^2). 

Assuming the source dust temperatures of 10 K, as in 
Section |4.2| are equal to the excitation temperature in the 
protostellar envelope, it follows from Equation |14| the ratio 
12 CO 6-5/ 12 CO 3-2 is -0.001 in the optically thin case. The 
ratio between the two lines is low due to the low temperature 
of the region, indicating there is less likelihood of detecting 
12 CO 6-5 in cooler regions of the cloud. The sources anal- 
ysed in NGC 1333, NGC 2071, and NGC 2024 with high 
12 CO 3-2 contamination correspond to molecular clouds at 
temperatures of 20 to 25 K. At 25 K, the ratio of 12 CO 6- 
5/ 12 CO 3-2 is 0.147 in the optically thin case. However, out- 
flows can contain even higher temperatures, ranging from 50 
to 150 K flHatchell et al.|1999| [van Kempen et al.|2009) . At 
50 K, the ratio is 0.769, indicating CO 6-5 is much more 
likely to be detected from outflows. If, on the other hand, 
both lines are optically thick, then the ratio tends to 1 as 



is known to be the case for 12 CO 3-2 jBuckle et al. 
Curtis et alT]|201Q ). 



2010 



Using the ratio 12 CO 6-5/ 12 CO 3-2, we can estimate 
typical peak fluxes for 12 CO 6-5: an excitation tempera- 
ture of 25 K and a source with a typical 12 CO 3-2 in- 
tegrated intensity of 100 K km s" 1 will produce a corre- 
sponding 12 CO 6-5 flux contribution of 8 mJy beam -1 for 
Grade 2 weather in the 450 um SCUBA-2 map. With the 
SCUBA-2 450 urn sensitivit y for the Gould Belt Survey 
( Ward-Thompson et al.|2007) at a RMS of -35 mJy beam" 1 
for Grade 2 weather, the CO 6-5 flux contribution would 
not be detected. Even in the case of optically thick emis- 
sion where the ratio 12 CO 6-5/ 12 CO 3-2 is 1, the 12 CO 6-5 
is estimated to be 57 mJy beam -1 , which is under the 5cr 
detection limit. 

Using published observations, 12 CO 6-5 contamination 
can be studied in further detail. 12 CO 6-5 data were taken 
for IR AS 2A, 4A, and 4B in NGC 1333 by |Yildiz et al.] 
(2010). These sources are particularly bright and possibly 
intermediate-mass protostars. 12 CO 6-5 integrated intensi- 
ties at the positions of the protostars were 57 K km s" 1 , 
122 K km s" 1 , and 43 K km s" 1 respectively. Using 
the 12 CO 6-5 conversion factors calculated in this study 
for Grade 2 weather, the corresponding CO contamina- 
tion to the SCUBA-2 450 urn dust continuum would be 
32 mJy beam" 1 (IRAS 2A), 70 mJy beam" 1 (IRAS 4A), 
and 25 mJy beam" 1 (IRAS 4B) for the CO contamination 
to the SCUBA-2 450 um dust continuum. Peak fluxes from 
SCUBA were 2355 mJy beam" 1 , 7000 mJy beam" 1 , and 
3025 mJy beam" 1 respectively jHatchell etaT]|2005| ) . The 
SCUBA-2 450 um dust emission peaks are a factor of several 
hundred times larger than the 12 CO 6-5 contribution. For 
these bright protostars, the CO contamination is insignifi- 
cant at 450 um. 



Dust continuum fluxes in the 450 um SCUBA-2 band 
are estimated to be a factor of 6 to 12 higher than fluxes in 
the 850 um band. In the Ray leigh- Jeans approximation, the 
dust optical depth increases as A _/3 with /3 between 1 and 2 
and the corresponding flux density increases as A" 3 to A" 4 
(Wilson 2009). The increase in continuum flux at 450 um 
clearly outweighs the expected contribution from 12 CO 6- 
5. At most the 12 CO 6-5 integrated main-beam emission 
is the same as the 12 CO 3-2 emission assuming optically 
thick emission. This situation corresponds to a molecular 
line ratio of 1 and a contamination flux ratio of 0.84 (as- 
suming Grade 2 weather). Since 450 um continuum fluxes 
increase with respect to the 850 um continuum, 450 um con- 
tinuum measurements would more likely be contaminated by 
strong molecular outflows instead of other means, i.e. nearby 
stars or ambient cloud emission. Nonetheless, potential con- 
tamination could occur in particularly low flux sources with 
nearby molecular outflows, such as the sources in NGC 2071 
that were near to confirmed outflow candidates (see Section 

Of. 

In the case of large-scale 12 CO 6-5 emission, it is also 
necessary to include the secondary beam in the calculation 
of the 450 um telescope beam (Equation [9}. The total beam 
area for 450 um increases by a factor of 2.0 (determined from 
an effective FWHM calculated in Section [Z2| . Using Equa- 
tion 11 with the the efficiency factor ?7f ss , the 12 CO 6-5 con- 
version factor would increase by a factor of 1.7. The increase 
in the expected 450 um dust continuum flux still exceeds 
any change in the 12 CO line conversion factors. Therefore, 
our conclusion that there is little CO contamination in the 
450 um maps, as discussed above, still holds. 



6 CONCLUSIONS 

In this study, the 12 CO line contamination factors for the 
450 um and 850 um SCUBA-2 continuum bands were cal- 
culated under different atmospheric conditions (weather 
grades 1 to 5). These contamination factors were then ap- 
plied to three different regions, NGC 1333, NGC 2071, and 
NGC 2024, in order to study the HARP 12 CO 3-2 flux con- 
tribution to the SCUBA-2 850 um measurements using a list 
of sources for each region. Sources with high 12 CO contam- 
ination (greater than 20 per cent) were analysed in further 
detail to determine the cause of the contamination. The fol- 
lowing can be concluded from this study: 

(i) For the 850 um SCUBA-2 filter profile, the 
12 CO 3-2 contamination factors increase as the sky 
opacity T225 increases. The contamination factors 
(mJy beam" 1 per K km s" 1 ) of 12 CO to the 850 um 
dust emission are, by weather grade, (Grade 1) 0.63; 
(Grade 2) 0.68; (Grade 3) 0.70; (Grade 4) 0.74; 
(Grade 5) 0.77. 

(ii) For the 450 um SCUBA-2 filter profile, the 
12 CO 6-5 contamination factors decrease as the sky 
opacity T225 increases due to the atmosphere transmis- 
sion steeply declining at higher opacity grades. The 
contribution factors (mJy beam" 1 per K km s" 1 ) of 
12 CO to the 450 um dust emission are, by weather 
grade, (Grade 1) 0.64; (Grade 2) 0.57; (Grade 3) 0.51; 
(Grade 4) 0.41; (Grade 5) 0.35. 
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(iii) The 12 C0 3-2 contribution to the 850 pun SCUBA-2 
dust continuum is typically under 20 per cent for all of the 
regions studied. However, in regions of molecular outflows, 
the 12 CO can reach a flux contribution of ~ 68 mJy beam -1 
for the sources studied, dominating the dust continuum 
in sources with both high and low continuum flux densi- 
ties (up to 500 mJy beam -1 ) with a contribution up to 
79 per cent contamination. Peak 12 CO fluxes in molecu- 
lar outflows in the regions reached even higher levels, up 
to 154 mJy beam -1 . There is strong evidence that 12 CO 3- 
2 contamination, while mostly minimal, is a major poten- 
tial source of confusion that can be observed directly in the 
850 \im dust continuum maps resembling protostellar cores 
or filamentary structure (as may be the case in NGC 1333, 
NGC 2071, and NGC 2024). 

(iv) Even though we have no 12 CO 6-5 molecular line 
maps to study in further detail, in hot (50 K) regions, 
e.g. molecular outflows, the ratio of main-beam tempera- 
ture T MB for 12 CO 6-5/ 12 CO 3-2 is -0.769. However, CO 
contamination to the 450 |^m source fluxes is not expected 
to be as much of an issue because of the expectation for the 
450 |^m dust emission to be a factor of 6 to 12 times brighter 
than the 850 um fluxes. 
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Table 3: List of sources used for the study of 12 CO contamination to the SCUBA-2 850 fim dust continuum. Source numbers correspond 
to the arbitrary number assigned to sources for this study. Source numbers corresponding to original studies (NGC 1333 sources obtained 
from |Hatchell et al.|2007| and NGC 2071 and NGC 2024 sources obtained from |Nutter fe Ward-Thompson|2007| are listed under 'Other 
Source ID'. 



Region 


Source 


RA 


Dec 


Other Source ID 






(J2000) 


(J2000) 




NGC 1333 


1 


03:29:10.4 


31 


13:30 


HRF41 




2 


03:29:12.0 


31 


13:10 


HRF42 




3 


03:29:03.2 


31 


15:59 


HRF43 




6 


03:28:55.3 


31 


14:36 


HRF44 




13 


03:29:01.4 


31 


20:29 


HRF45 




14 


03:29:11.0 


31 


18:27 


HRF46 




17 


03:28:59.7 


31 


21:34 


HRF47 




18 


03:29:13.6 


31 


13:55 


HRF48 




20 


03:28:36.7 


31 


13:30 


HRF49 




21 


03:29:06.5 


31 


15:39 


HRF50 




22 


03:29:08.8 


31 


15:18 


HRF51 




24 


03:29:03.7 


31 


14:53 


HRF52 




26 


03:29:04.5 


31 


20:59 


HRF53 




29 


03:29:10.7 


31 


21:45 


HRF54 




30 


03:28:40.4 


31 


17:51 


HRF55 




32 


03:29:07.7 


31 


21:57 


HRF56 




33 


03:29:18.2 


31 


25:11 


HRF57 




35 


03:29:16.5 


31 


12:35 


HRF59 




38 


03:28:39.4 


31 


18:27 


HRF60 




39 


03:29:17.3 


31 


27:50 


HRF61 




42 


03:29:07.1 


31 


17:24 


HRF62 




43 


03:29:18.8 


31 


23:17 


HRF63 




46 


03:29:25.5 


31 


28:18 


HRF64 




47 


03:29:00.4 


31 


12:02 


HRF65 




53 


03:29:05.3 


31 


22:11 


HRF66 




55 


03:29:19.7 


31 


23:56 


HRF67 




62 


03:28:56.2 


31 


19:13 


HRF68 




63 


03:28:34.4 


31 


06:59 


HRF69 




66 


03:29:15.3 


31 


20:31 


HRF70 




73 


03:28:38.7 


31 


05:57 


HRF71 




77 


03:29:19.1 


31 


11:38 


HRF72 




81 


03:28:32.5 


31 


11:08 


HRF74 




83 


03:28:42.6 


31 


06:10 


HRF75 




89 


03:29:04.9 


31 


18:41 


Bolo44 




94 


03:28:32.7 


31 


04:56 


Bolo26 


NGC 2071 


51 


05:47:23.7 


00 


11:02 


BN-547237+01102 




52 


05:47:06.8 


00 


12:30 


BN-547068+01230 




53 


05:47:10.6 


00 


13:18 


BN-547106+01318 




54 


05:47:05.1 


00 


13:21 


BN-547051+01321 




55 


05:47:05.0 


00 


14:49 


BN-547050+01449 




56 


05:47:23.9 


00 


15:07 


BN-547239+01507 




( 


n K . A T. 1 1 A 


00 


15:37 


Jd1N-o4 ( lz4+Ul0o ( 




58 


05:47:10.4 


00 


15:53 


BN-547104+01553 




59 


05:47:19.9 


00 


16:03 


BN-547199+01603 




60 


05:47:04.8 


00 


17:07 


BN-547048+01707 




61 


05:47:01.5 


00 


17:55 


BN-547015+01755 




62 


05:47:08.7 


00 


18:17 


BN-547087+01817 




63 


05:47:15.2 


00 


18:30 


BN-547152+01830 




64 


05:47:25.3 


00 


18:48 


BN-547253+01848 




65 


05:47:33.6 


00 


19:02 


BN-547336+01902 




66 


05:46:28.3 


00 


19:28 


BN-546283+01928 




67 


05:47:03.4 


00 


19:50 


BN-547034+01950 




68 


05:47:26.7 


00 


19:53 


BN-547267+01953 




69 


05:47:37.7 


00 


20:01 


BN-547377+02001 




70 


05:46:57.6 


00 


20:09 


BN-546576+02009 
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Table 3: continued from previous page. 



Region 


Source 


RA 
(J2000) 


Dec 
(J2000) 


Other Source ID 




71 


05:46:29.4 


00:20:10 


BN-546294+02010 




72 


05:47:34.9 


00:20:20 


BN-547349+02020 




73 


05:47:32.5 


00:20:26 


BN-547325+02026 




74 


05:47:01.0 


00:20:42 


BN-547010+02042 




75 


05:47:25.2 


00:20:59 


BN-547252+02059 




76 


05:47:10.3 


00:21:12 


BN-547103+02112 




77 


05:46:28.7 


00:21:14 


BN-546287+02114 




78 


05:47:16.0 


00:21:23 


BN-547160+02123 




79 


05:47:04.1 


00:21:58 


BN-547041+02158 




80 


05:46:25.3 


00:22:20 


BN-546253+02220 




81 


05:46:52.8 


00:22:23 


BN-546528+02223 




82 


05:47:11.9 


00:22:23 


BN-547119+02223 




83 


05:47:19.7 


00:22:31 


BN-547197+02231 




84 


05:47:06.9 


00:22:39 


BN-547069+02239 




85 


05:47:17.5 


00:22:40 


BN-547175+02240 




86 


05:46:59.1 


00:22:59 


BN-546591+02259 




87 


05:47:12.4 


00:23:11 


BN-547124+02311 




88 


05:47:06.7 


00:23:14 


BN-547067+02314 




89 


05:46:54.7 


00:23:24 


BN-546547+02324 




90 


05:47:10.4 


00:23:27 


BN-547104+02327 




91 


05:47:08.9 


00:23:56 


BN-547089+02356 




92 


05:46:57.2 


00:23:56 


BN-546572+02356 




93 


05:46:34.7 


00:23:59 


BN-546347+02359 




94 


05:46:58.0 


00:24:26 


BN-546580+02426 




95 


05:47:01.7 


00:24:52 


BN-547017+02452 




96 


05:46:25.7 


00:24:56 


BN-546257+02456 




97 


05:47:08.0 


00:25:05 


BN-547080+02505 




98 


05:46:45.9 


00:25:07 


BN-546459+02507 




99 


05:47:01.4 


00:26:14 


BN-547014+02614 




100 


05:46:38.0 


00:26:53 


BN-546380+02653 


NGC 2024 


33 


05 42 03.0 


-02 04 23 


BS-542030-20423 




34 


05 42 10.3 


-02 04 20 


BS-542103-20420 




35 


05 42 03.5 


-02 02 24 


BS-542035-20224 




36 


05 41 57.1 


-02 01 00 


BS-541571-20100 




37 


05 41 52.9 


-02 00 21 


BS-541529-20021 




38 


05 41 49.3 


-01 59 38 


BS-541493-15938 




40 


05 42 00.0 


-01 58 01 


BS-542000-15801 




41 


05 41 49.1 


-01 58 03 


BS-541491-15803 




42 


05 41 45.2 


-01 56 31 


BS-541452-15631 




43 


05 41 35.4 


-01 56 29 


BS-541354-15629 




45 


05 41 44.5 


-01 55 39 


BS-541445-15539 




46 


05 41 44.2 


-01 54 43 


BS-541442-15443 




47 


05 41 32.1 


-01 54 26 


BS-541321-15426 




48 


05 41 19.9 


-01 54 16 


BS-541199-15416 




49 


05 41 42.0 


-01 53 59 


BS-541420-15359 




51 


05 41 44.2 


-01 52 41 


BS-541442-15241 




52 


05 41 16.6 


-01 51 19 


BS-541166-15119 




53 


05 41 36.7 


-01 51 06 


BS-541367-15106 




54 


05 41 23.4 


-01 50 27 


BS-541234-15027 




55 


05 41 32.9 


-01 49 53 


BS-541329-14953 




56 


05 41 36.4 


-01 49 24 


BS-541364-14924 




57 


05 41 27.6 


-01 48 13 


BS-541276-14813 




58 


05 41 11.3 


-01 48 12 


BS-541113-14812 




59 


05 41 13.3 


-01 47 35 


BS-541133-14735 







20 E. Drabek et al 

APPENDIX A: MOLECULAR OUTFLOW 
SPECTRA CRITERIA 

The full version of this Appendix is available as Supporting 
Information to the online version of this article. 
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Figure Al. 12 CO 3-2 spectra for sources further analysed in Section [o] regarding the molecular outflow analysis. The core velocities are 
listed in Section liXTl for NGC 1333 and in Table[|]for NGC 2071 and NGC 2024. Parameters ±3 km s" 1 for NGC 1333 and ±4 km s" 1 
for NGC 2071 and NGC 2024 at 1.5 K were used to classify the presence of outflows. 



